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ABSTRACT 

We present interferometric observations with resolution of ~ 3" of the iso- 
lated, low-mass protostellar double cores CG30 and BHR71 in the N2H+ (1 — 0) 
line and at 3 mm dust continuum, using the Australian Telescope Compact Array 
(ATCA). The results are complemented by infrared data from the Spitzer Space 
Telescope. In CG30, the 3 mm dust continuum images resolve two compact 
sources with a separation of ~21'/7 (~8700AU). In BHR71, strong dust contin- 
uum emission is detected at the position of the mid- infrared source IRSl, while 
only weak emission is detected from the secondary mid-infrared source IRS2. 
Assuming standard gas to dust ratio and optically thin 3 mm dust continuum 
emission, we derive hydrogen gas masses of 0.05 — 2.1 Mq for the four sub-cores. 
N2H''" (1 — 0) line emission is detected in both CG 30 and BHR71, and is spatially 
associated with the thermal dust continuum emission. By simultaneously fitting 
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the seven hyperfine line components of N2H"'", we derive the velocity fields and 
find symmetric velocity gradients in both sources. Assuming that these gradients 
are due to core rotation, we estimate the specific angular momenta and ratios of 
rotational energy to gravitational energy for all cores. Estimated virial masses of 
the sub-cores range from 0.1 — 0.6 Mq. We also find that the N2H"'" emission is 
strongly affected by the outflows, both in terms of entrainment and molecule de- 
struction. Spitzer images show the mid-infrared emission from all four sub-cores, 
which is spatially associated with the 3 mm dust continuum emission. All four 
sources appear to drive their own outflows, as seen in the shock-excited 4.5 /im 
images. Based on the ATCA and Spitzer observations, we construct spectral en- 
ergy distributions (SEDs) and derive temperatures and luminosities for all cores. 
The analysis of the SEDs suggests an evolutionary discrepancy between the two 
sub-cores in both CG30 and BHR71, which could be due to effects of relative in- 
clinations. Based on the morphology and velocity structure, we suggest that the 
sub-cores in CG 30 were formed by initial fragmentation of a fllamentary prestel- 
lar core, while those in BHR 71 could originate from rotational fragmentation of 
a single collapsing protostellar core. 

Subject headings: binaries: general — ISM: globules — ISM: individual (CG 30 
and BHR 71) — ISM: kinematics and dynamics — stars: formation 

1. INTRODUCTION 

Although statistical properties of binary stars have been determined over the past two 
decades (see e.g., Reipurth et al. 2007 for recent reviews), many key questions concerning 
their origin are still poorly understood. What is the formation mechanism of binary /multiple 
systems? How are mass and angular momentum distributed during their formation? What is 
the difference between cores forming binaries and those forming single stars? To answer these 
questions, direct observations of the earliest, deeply embedded phase of binary star formation 
are needed. This phase is unfortunately not accessible to optical and near-infrared (NIR) 
wavelengths, due to the large amounts of circumstellar material present. Observations of the 
gas and optically thin dust emission at millimeter (mm) wavelengths are therefore needed, to 
probe the system kinematics and individual envelope masses. However, these observations 
were long hampered by the low angular resolution of mm telescopes, and only recently have 
the earliest phases of binary star formation been observationally identifled and studied in 
detail thanks to the availability of large (sub-) mm interferometers, although the number of 
known systems is still very small (Looney et al. 2000; Launhardt 2004). 
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To search for binary protostars and to derive their kinematic properties, we have started 
a systematic program to observe, at high angular resolution, a number of isolated low-mass 
prcstellar and protostellar molecular cloud cores. The initial survey was conducted at the 
Owens Valley Radio Observatory (OVRO) mm array (Launhardt 2004; Chen, Launhardt, & 
Henning 2007, hereafter Paper 1; Launhardt et al., in prep.), and is now continued with the 
Australia Telescope Compact Array (ATCA) and the IRAM Plateau de Bure Interferometer 
(PdBI) array. In this paper we present ATCA observations of two southern protobinaries in 
the N2H"'" (1 — 0) molecular line and at 3 mm dust continuum, together with complementary 
mid-infrared (MIR) data from the Spitzer Space Telescope (hereafter Spitzer). 

CG30 (also known as BHR12 or DC 253.3—1.6) is a bright-rimmed cometary globule 
located in the Gum Nebula region. The distance towards CG30 is somewhat uncertain, 
with estimates ranging from 200 pc (Knude et al. 1999) to 400 pc (Brandt 1971; Reipurth 
1983). For consistency with earher papers (e.g., Henning & Launhardt 1998), we use here 
400 pc. The globule harbors an elongated protostellar core as seen in the single-dish mm 
dust continuum image (Launhardt et al. in prep.). Higher resolution submm continuum 
observations (SCUBA) resolve the source into two sub-cores with a projected separation of 
~20" (~8000AU) and masses of 0.17 ±0.05 and 0.14 ±0.05 Mq (Henning et al. 2001). The 
northern core is associated with a NIR source, which drives the Herbig-Haro flow HH 120 
(see Hodapp & Ladd 1995 and references therein). The newly discovered southern core is 
the origin of a protostellar jet with position angle (P. A.) 44° (Hodapp & Ladd 1995), but no 
NIR source is seen at this position (see Launhardt et al. 2001). 

BHR71 (also known as DC 297.7—2.8) is an isolated Bok globule located at a distance 
of ~ 200 pc (Bourke et al. 1997; hereafter B97). A highly-coUimated bipolar outflow, which 
is lying almost in the plane of the sky, was discovered by CO observations in this region. The 
driving source is associated with IRAS 11590—6452 and was classified as a Class protostar 
with a total luminosity of ~ 9L0 (B97). ISOCAM observations have revealed that the 
IRAS source is associated with two embedded protostars, IRSl and 1RS2, with a projected 
separation of ~ 17" (~3400AU; Bourke 2001; hereafter BOl). IRSl and 1RS2 each drive 
a CO outflow: the well-known large-scale collimated bipolar outflow is driven by IRSl and 
another fainter and smaller bipolar outflow is driven by IRS2 (see BOl and Parise et al. 2006). 
Only IRSl appears to be associated with a substantial amount of circumstellar material, but 
neither is directly detected at NIR wavelengths (BOl). 
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2. OBSERVATIONS AND DATA REDUCTION 

2.1. ATCA Observations 

Millimeter interferometric observations at 95 GHz of CG 30 and BHR 71 were carried out 
using ATCA with five 22 m telescopes in May and August 2005. Observations were obtained 
in two different array configurations (H168 and H75) with projected baselines ranging from 
22 to 180 m. All antennas were equipped with cooled SIS receivers, which provided average 
system temperatures of 200 — 350 K at the observing frequency. A digital correlator was 
used with 2 independent spectral windows. The narrow window (bandwidth ~8MHz), with 
a channel width of 0.019 MHz, was centered on the N2H+ (1-0) line at 93.17GH^ The 
broad window (bandwidth ~ 128 MHz) was centered at 95.0 GHz and was used to measure 
the 3.1mm dust continuum emission simultaneously with N2H"'". The two sources were 
observed with 2-point mosaics each. The primary beam size at 93 GHz is ~ 38". Amplitude 
and phase were calibrated through frequent observations of quasars nearby to each source 
(0745—330 for CG30 and 1057—797 for BHR 71), typically every 20 minutes, resulting in an 
absolute position uncertainty of < 0'.'2. Flux densities were calibrated using the secondary 
calibrator 1253—055, the flux of which was regularly compared to Uranus and adopted as 
19.0 Jy for May observations (H168 configuration) and 14.7 Jy for August observations (H75 
configuration). Additional effort was made to improve the gain-elevation calibration of the 
antennas, which can significantly affect the flux density scale, especially when observing at 
high elevation. The estimated total flux uncertainty is < 20%. Observing parameters are 
summarized in Table 1. 

The data were calibrated and images produced using MIRIAD (Sault et al. 1995) 
and its CLEAN algorithm, with "robust" uv weighting parameter +1 (Briggs et al. 1999). 
Synthesized beam sizes are 3" — 4". Noise levels (la rms) in the final maps are 0.5 — 2 
mJy/beam for the continuum and 20 — 65 mJy/beam for the N2H"'" line (see Table 1). Further 
analysis and figures were done with the GILDAS[^ software package. 



^During the observations the sky frequency changed by less than one channel due to the lack of doppler 
tracking at ATCA, and corrections were applied offline to obtain correct frequencies and LSR velocities 

2http://www.iram.fr/IRAMFR/GILDAS 
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2.2. Spitzer Observations 

Mid-infrared data of CG 30 and BHR 71 were obtained from the Spitzer Science Centeij^ 
CG30 was observed on 2004 April 9 with the Multiband Imaging Photometer for Spitzer 
(MIPS: AOR key 9426688) and May 26 with the Infrared Array Camera (IRAC; AOR key 
5097216). BHR 71 was observed on 2004 June 10 with IRAC (AOR key 5107200) and 2005 
March 7 with MIPS (AOR key 9434112). Both sources were observed as part of the c2d 
Legacy program (Evans et al. 2003). 

The data were processed by the Spitzer Science Center using their standard pipeline 
(version S14.0) to produce Post Basic Calibrated Data (P-BCD) images, which are fiux- 
calibrated into physical units (MJy sr~^). Flux densities in the IRAC bands were measured 
with aperture photometry in the IRAF APPHOT package, using the radii, background 
aperture annuli, and aperture corrections recommended by the Spitzer Science Center. The 
results were compared to c2d, which used PSF fitting, and found to be within the uncertain- 
ties. Flux densities in the MIPS bands were measured with GILDAS because sources in the 
MIPS images are not fully resolved (see §3.3). Further analysis and figures were completed 
with GILDAS. 

3. RESULTS 

3.1. Dust Continuum 

The 3 mm dust continuum image of CG 30 (Fig. la) shows two compact sources with an 
angular separation of 21'.'7 ± 0'.'6, corresponding to a projected linear separation of 8700 ± 240 AU 
at a distance of 400 pc. Following Henning et al. (2001), we refer to the northern source as 
CG30N and to the southern source as CG30S. From Gaussian uv plane fitting, we derive 
flux densities of 15.8 ± 3.2 m J}|^ for source N and 6.0±1.3mJy for source S. The large-scale 
common envelope, detected in the submm single-dish maps with a radius of ~ 14000 AU and 
a flux density of ~7.4Jy (Henning et al. 2001), is resolved out by the interferometer at 
3 mm. Source positions and deconvolved FWHM sizes of the two embedded sources, derived 
from Gaussian uv plane fitting, are listed in Table 2. 

In the 3 mm dust continuum image of BHR 71 (Fig. lb), strong emission is detected at 



■^http:/ /ssc. spitzer. caltech.edu 

^The error bar is derived from V c^aii ~^ '^Rt ' where acaii is the uncertainty from cahbration (~ 20% of flux 
density) and afn is the uncertainty from Gaussian fitting. 
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the position of IRSl, and only weak emission (~3(7 level) is detected at the position of 
IRS2. The flux densities of IRSl and IRS2 are derived to be 140 ± 28 mJy and 2.8 ± 2.1 mJy, 
respectively. The large-scale envelope detected in the 1.3 mm single-dish map, with a radius 
of ~9000AU and a flux density of ~3.7Jy (B97), is also resolved out here. Positions 
and FHWM sizes of the sources are listed in Table 2. The angular separation of 17" ± 1" 
between IRSl and IRS2 corresponds to a projected linear separation of 3400 ± 200 AU at a 
distance of 200 pc. We also note that IRSl is elongated northwest-southeast and consists of 
two separate peaks in the region enclosed by the 5(7 level (see Fig. lb). The main peak is 
spatially coincident with the MIR source and the fainter peak is located ~ 2" southeast of 
IRSl (see below §3.3). 

Assuming that the 3 mm dust continuum emission is optically thin, the hydrogen gas 
mass Mh = M(H) -|- 2 M(H2) in the circumstellar envelope (excluding Helium) was calcu- 
lated with the same method as described in Launhardt & Henning (1997). We adopt an 
interstellar hydrogen-to-dust mass ratio of 110, and a dust opacity 0.2cm^g ^ (us- 

ing Ki 3inin = 0.8cm^g~^ and k. oc u^'^), a fairly typical value for dense protostellar cores 
(Ossenkopf & Henning 1994). Dust temperatures are derived from SED fltting (§4.1) and 
are listed in Table 7. The derived hydrogen gas masses (0.05 — 2.1 Mq), together with mean 
volume densities (0.4 — 2.6 x 10''cm~'^) and column densities (1.5 — 9.9 x 10^'^ cm~^), are 
hsted in Table 2. The resulting optical depths are Tsmm ~ 0.4 — 3 x 10~^, thus justifying the 
optically thin approximation. 



3.2. N2H+(l-0) 

N2H"'' emission is detected from both CG30 and BHR71. Figure 2a shows the velocity- 
integrated N2H+ intensity image of CG30. Two cores, spatially associated with the 3 mm 
dust continuum sources, are seen. The northern core is elongated east- west with a long 
(~ 20") extension to the west, along the direction of the protostellar jet HH 120 (see Fig. 7c). 
The southern core is more compact and peaks at the position of the dust continuum source. 
The mean radii of the N2H''" cores (Table 4) were measured with the same method described 
in Paper I. A larger-scale N2H"'" cloud core, detected in the Mopra single-dish map with a 
radius of ~ 8000 AU (P. Barnes et al. in prep.), peaks between the two sub-cores (see Fig. 2a) 
and is resolved out by the interferometer (more than 90% flux is missing). 

Figure 2b shows the integrated N2H"'" intensity image of BHR 71. Two cores are found to 
the east and west of IRSl (see below and discussion in § 4.4). We refer to these as BHR71E 
and BHR 71W, respectively. The two cores are elongated in the north-south direction. Sev- 
eral smaller clumps are also seen north and south of the two main cores, along both sides 
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of the large-scale CO outflow (see Fig. 2b). For BHR71, the Mopra N2H"'" map (P. Barnes 
et al. in prep.) again shows a large-scale cloud core with one peak (~ 15" offset IRSl), and 
does not line up well with the two dust continuum sources (see Fig. 2b). 

Figure 3 shows the N2H"'" spectra at the peak positions of CG 30 and BHR 7l|^ The 
spectra were fitted using the hyperfine program in CLASS. The fitting results, such as LSR 
velocities (Vlsr), intrinsic hne width (Af ; corrected for instrumental effects), total optical 
depths (rtot), and excitation temperatures (Tex), are listed in Table 3. 

Figure 4 shows the mean velocity fields of CG 30 and BHR 71, derived from the N2H"'" line 
maps with the fitting routine described in Paper I. The jet/outfiow information is also shown 
in each map. In CG 30, the southern core shows a well-ordered velocity field, with gradient 
parallel to the outflow direction. The northern core shows a more complicated velocity field, 
but the gradient in the inner core is also parallel to the outflow direction. In BHR 71, there 
seems to be a general velocity gradient across the two N2H"*' cores, which is approximately 
perpendicular to the axis of the large-scale CO outflow. This may indicate that the two cores 
are actually part of one physical structure associated with IRSl (see discussion in §4.4). A 
least-squares fitting of the velocity gradients has been performed using the routine described 
in Goodman et al. (1993). The results are summarized in Table 5 and discussed in §4.2. 

Figure 5 shows the spatial distribution of N2H"'" line widths for both sources. The line 
widths are roughly constant within the interiors of the cores, which is consistent with the 
observational results in Paper I. The mean line widths were derived through Gaussian fitting 
to the distribution of line widths versus solid angle area in the maps (see Fig. 6). We find 
that the sub-cores in each object have roughly equal line width, but the mean line width in 
CG30 (~0.5kms-^) is ~ 1.7 times larger than that in BHR 71 (~0.3kms-^). 

Assuming that the observed N2H"'" line widths are not dominated by systematic gas 
motions, the virial mass of the cores has been calculated as: 

5 RAvl,^ 
81n2 a^irG 

where G is the gravitational constant, R is the FWHM core radius, and Afave is the line 
width of the emission from an "average" particle with mass m^^e = 2.33 amu (assuming gas 
with 90% H2 and 10% He). The coefficient Ovir = (1 ^ p/3)/(l — 2^/5), where p is the 
power-law index of the density profile, is a correction for deviations from constant density 



^During the observations towards BHR 71, the correlator was not well centered due to an uncertainty in 
the Doppler correction calculation, resulting in the N2H+ JFiF ^ 101 — 012 component not being covered. 
However, the absence of this line component did not affect our final results. 
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(see Williams et al. 1994). In our calculations, we assume p — 1.5 (see Andre et al. 2000) 
and use CKyir = 1.25. Avgye is derived from the observed spectra by 

A^L = + 81n2^(-i- - -J-), (2) 

mn mave "7.obs 

where Afobs is the observed mean line width of N2H+ and mobs is the mass of the emitting 
molecule (here we use mN2H+ = 29 amu). We derive virial masses between 0.1 and 0.6 M©. 
The results are listed in Table 4. 

The N2H+ column density has been calculated independently from the line intensity 
using the equation given by Benson et al. (1998): 

7V(N2H+) = 3.3 X 10" ^^^^g?^ (cm-^), (3) 

where r is the total optical depth, Av is the intrinsic line width in kms~^, and Tex is the 
excitation temperature in K. The gas-phase N2H"'" mass of the core was then calculated from 
^^N2H+ ~ ^(N2H+)pcak X '"^N2H+ X c?^ X l^pwHM, where d is the distance from the Sun and 
^FWHM is the solid angle enclosed by the FWHM contours for each core. 

Assuming that the gas mass and virial mass derived from the N2H"'" data are the same, 
we derived the average fractional abundance of N2H+ in each core (sec Table 4). The average 
value (X(N2H+)) ~3.0 x IQ-io for CG30 and BHR71 is close to the mean value found in 
Paper I (~ 3.3 x 10"^") for nine other protostellar cores. 



3.3. Spitzer Images 

Figure 7 shows the Spitzer images of CG30. The infrared emission from CG30N and 
CG30S is detected at all IRAC bands (3.6//m — 8.0 /xm). Fig. 7a shows a wide-field IRAC 
band 2 (4.5 ;um) image. Centered at GG30S is a highly collimated bipolar jet, with P.A. 
~40°. The knots in the jet are labeled with the same numbers as in Hodapp & Ladd (1995). 
The most distant knot (No. 8) is ~ 90" away from CG 30S. Assuming a typical jet speed of 
100 km s~^ (Reipurth & Bally 2001), an inclination angle of 90°, and a distance of 400 pc, 
the dynamical age of the jet is estimated to be ~ 1700 yr. CG 30N appears to be the driving 
source of HH 120, which is ~ 5" in size and extends to the west. Knot No. 6, located to the 
east of CG 30N, is probably ejected by CG 30N and part of the same outflow as HH 120. 

Figs. 7b and 7c show enlarged views of the two sources, overlaid with the contours from 
the ATCA 3 mm dust continuum and N2H+ images. The two infrared sources are spatially 
coincident with the 3 mm dust continuum and N2H+ sources. However, when viewed in detail. 



-9- 



CG 30S is elongated at the infrared bands and the continuum source is located at the apex 
of the infrared emission, implying that the infrared emission at IRAC bands from CG 30S is 
due to scattered light in a cavity evacuated by the jet/outflow. In contrast, CG SON shows a 
point-like structure at all IRAC bands coincident with the circumstellar mm dust emission 
peak, suggesting that the source is directly detected at NIR wavelengths (A < 5/im). The 
N2H+ emission from CG SON spatially follows the direction of the protostellar jet and the 
long extension to the west matches exactly the HH 120 flow (see Fig. 7c), indicating that the 
jet has a strong effect on the morphology of the N2H+ emission. 

In the Spitzer MIPS 1 (24 /im) image shown in Fig. 7d, CG SO is again resolved in two 
sources, but the emission is dominated by CG SON and only weak emission is found at the 
position of CG30S. In the MIPS 2 (70 /xm) image (see Fig. 7e), the two sources are not fully 
resolved, but two peaks, with flux ratio ~ 2:1, can be clearly distinguished. Flux densities 
of CG SON and CG SOS in the IRAC and MIPS bands are measured (see § 2.2) and hsted in 
Table 6. 

The Spitzer images of BHR 71 are shown in Fig. 8, with the same sequence as in Fig. 7. 
The infrared emission from IRSl and IRS2 is detected at all IRAC bands. A large-scale (~ 
160" in length) bipolar jet, centered at IRSl with a P. A. of 165°, is seen in the IRAC images 
(Fig. 8a). The northern jet, spatially coincident with the red-shifted CO outflow, is S-shaped, 
while the southern jet, containing the HH object HHS21 (Corporon & Reipurth 1997), shows 
a V-shaped structure at the apex. This V-shaped structure may represent a conical cavity 
evacuated by the successive bow-shocks traced by the infrared emission (Fig. 8a) and the 
blue-shifted CO outflow (see B97 and Parise et al. 2006). Another bipolar jet, at P. A. ~ 
S0°, is found with IRS2 being in the center. Its northwest lobe, containing another HH 
object HHS20 (Corporon & Reipurth 1997), also shows a V-shaped structure at the apex 
and could be explained in the same way. 

IRSl and 1RS2 are spatially coincident with the dust continuum sources detected with 
ATCA (see Fig. 8b). We note that the IRS2 dust continuum source is located at the apex 
of the infrared emission and could be explained in the same way as CG SOS. The elongated 
structure and secondary peak found in the ATCA dust continuum image match the left wall 
of the outflow cavity, suggesting they result from the jet/outflow action (for a similar case, 
see Gueth et al. 200S). The N2H"'" emission is located on both sides of the large-scale CO 
outflow and basically matches the wall of the cavity (see Fig. 8c). At the MIPS 1 band, 
BHR 71 is barely resolved into two sources and the emission is dominated by IRSl (sec 
Fig. 8d). The MIPS 2 image does not resolve the two sources and the emission is peaked at 
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the position of IRSlj^ (see Fig. 8e). Flux densities of IRSl and IRS2 are listed in Table 6. 

4. DISCUSSION 

4.1. Spectral Energy Distributions and Evolutionary Stages 

Figure 9 shows the spectral energy distributions (SEDs) of CG30 N and S and BHR71 
IRSl and IRS2, based on the infrared (ISOCAM, Spitzer, and IRAS), sub-mm (SCUBA, 
available only for CG30), and mm (SEST and ATCA) observations. The NIR data of 
CG30N are adopted from Persi et al. (1990). The SCUBA and SEST data for CG30 
are adopted from Henning et al. (2001) and Henning & Launhardt (1998), respectively. 
The ISOCAM and SEST data for BHR71 are adopted from BOl and B97, respectively 
Here we do not explicitly list all flux values, but show graphically the SEDs. Since IRAS 
observations could resolve neither CG30 nor BHR71, flux ratios at the IRAS wavelengths of 
10:1 (CG30N:CG30S) and 20:1 (BHR71 IRSl : IRS2) were inferred from the Spitzer and 
ATCA observations. 

In order to derive luminosities and bolometric temperatures, we first interpolated and 
then integrated the SEDs, always assuming spherical symmetry. Interpolation between the 
flux densities was done by a grey-body fit to all points at A > 100 /iuQ using 

= B,{T^){1 - e-^-')n, (4) 

where ByiT^) is the Planck function at frequency p and dust temperature Tj, Ty is the dust 
optical depth as a function of frequency r cx z/^-^, and Vt is the solid angle of the source. 
A simple logarithmic interpolation was used between all points at A < 100 /xm. The fitting 
results, such as dust and bolometric temperatures, sub-mm (A > 350 /xm) and bolometric 
luminosities, are listed in Table 7. 

Based on these results, we try to address the evolutionary stages of CG30 and BIIR71. 
A detailed definition and discussion for early stellar evolutionary phases can be found in 
Andre et al. (2000) and Froebrich (2005). The i^submm/-Z^boi ratios of all sources are ^ 
0.5% (the standard boundary of Class protostars, see Andre et al. 2000)r] and the four 



^Thc offset (^3") between the MfPS2 emission peak; and the 3 mm emission peak; is much smaller than 
the FWHM of MfPS2 PSF (> 10"), and is not significant. 

^The 3 mm points were ignored in the fitting to CG 30 and BHR 71 IRS2 to give higher priority to the 
submm data, resulting in much better fitting. 

^Our LsubmrnZ-^boi ratios are larger than those found by Froebrich (2005). We attribute this to the fact 
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sources each drive a bipolar jet (see § 3.3). However, the bolometric temperature of CG30N 
is ~ 100 K and the object is also directly detected at NIR wavelengths, suggesting CG SON 
is a Class I young stellar object. In contrast, the low bolometric temperature (37 K) of 
CG30S suggests it is a Class protostar. In BIIR71, both IRSl and IRS2 have bolometric 
temperatures less than 70 K (see Table 7). Nevertheless, IRSl might be directly detected at 
NIR wavelengths (see Fig. 8a), suggesting that it is a transition object between Class and 
I, while IRS2 could be a Class protostar. 

It must be noted that the analysis above does not take into account inclination effects: 
considering a protostar embedded in a circumstcllar disk/envelope, its infrared emission 
could be detected through the outflow cavity when the system is face-on, but is not seen 
when it is edge-on. In BHR71, the bipolar CO outflow powered by IRSl is lying roughly in 
the plane of sky, implying the latter case; the bipolar outflow driven by IRS2 appears to favor 
the same situation (see Parise et al. 2006). In CG 30, however, the relative inchnations could 
not be easily distinguished because the information about molecular outflows is still missing. 
The SED-based classification discussed above thus might be a result of both evolutionary 
stage and inclination. In particular in CG 30 we cannot disentangle the two effects. It is well 
possible that the different SEDs (and bolometric temperatures) reflect actually inclination 
effects rather than a difference in evolutionary stage. 



4.2. Gas Kinematics 

The thermal contribution to the NoH"*" line width is calculated by Au^^ = 81n2-^^, 

^ th "lobs' 

where k is the Boltzmann constant, Tk is the kinetic gas temperature, and mobs is the mass 
of the observed molecule. Assuming that at the high densities of > 10^ cm~^ (see Table 2) the 
kinetic gas temperature is equal to the dust temperature derived in §4.1 (~20K), the non- 
thermal contributions to the line widths (Ai;nt = ^/Au^gj^j^ — Av^^) were then calculated 
to be ~0.5kms-^ in CG30 and ~0.2kms-^ in BHR71 (see Table 3). These non-thermal 
line widths suggest that turbulence, the main contribution to the non-thermal line width 
(Goodman et al. 1998), cannot be ignored in the protostellar cores. On the other hand, 
the thermal line width of an "average" particle of mass 2.33 mn (assuming gas with 90% II2 
and 10% He), which represents the local sound speed, is ~0.62kms~^ at 20 K. The derived 
non-thermal contributions to the N2H"'" line width in both CG 30 and BHR 71 are smaller 



that Frobrich (2005) assumed the two objects were single cores, but we resolved them as binaries. Further- 
more, we have more data points at submm wavelengths (for CG 30) and high-resolution interferometric data 
points at 3 mm, which were all not available to Proebrich (2005). 
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than this local sound speed (i.e., the turbulent motion is subsonic). We also note that the 
mean line widths derived for BHR71 (~0.3kms~^) arc three times smaller than measured 
by single-dish observations in Mardones et al. (1997; ~0.9kms~^). Taking into account the 
systematic velocity variation across the core (~0.3kms~^; see Table 5), the combined line 
width in our maps is still smaller than the result from single-dish observations. It means that 
high-level (supersonic) turbulence occurs mainly in the extended envelope which is resolved 
out by the interferometer, but the inner core is much more "quiescent" . This is consistent 
with what we found in Paper I, namely that non-thermal motions are quickly damped from 
large-scale to smaller inner cores (see e.g.. Puller & Myers 1992). 

The velocity fields of CG30N, CG30S, and BHR71 show systematic velocity gradients 
(see Fig. 4). As discussed in Paper I, systematic velocity gradients are usually dominated 
by either rotation or outflow. In CG30, the gradients in both cores are parallel to the 
jets. Although there is no molecular outflow information available yet for CG30, these 
gradients are likely the results of outflows and we treat them as upper limits of underlying 
rotation velocity gradients. In BHR71, the velocity gradient across the two N2H"'" cores 
is roughly perpendicular to the axis of the large-scale CO outflow and could be explained 
by rotation. [Here we assume that the two cores are associated with IRSl (see §4.4).] 
The velocity gradients measured in CG30N, CG30S, and BHR71 are < 24.4, < 17.5, and 
7.8±0.5kms~^pc~^, respectively (see Table 5). Both the velocity gradient in BHR71 as 
well as the upper hmits for CG 30N and CG 30S are consistent with those found in Paper I. 

Assuming that the velocity gradients summarized in Table 5 are due to core rotation, the 
specific angular momentum J /M of the objects was calculated with the following equation: 

J/M = a,^,uR' = llp^^.R' ^ \ gR\ (5) 
6b — p Sim 7 

where the coefficient ctrot = | |E|) P is the power-law index of the radial density profile (p 
— 1.5; see §3.2), g is the velocity gradient, and i is the inclination angle to the line of sight 
direction (here we assume sini = 1). The derived J/M for CG30N, CG30S, and BHR71 
are listed in Table 5. It should be noted that for CG30 we derive only upper limits. The 
ratio of rotational energy to the gravitational potential energy was calculated by /9rot = 

^ 0.19 where Erot = \ ^ = \ "rot MR^ uJ^ and Egrav = | avuGM'^/R (the masses 
and radii used in the equations are virial masses and radii listed in Table 4) . The estimated 
/?rot values for CG30N, CG30S, and BHR71 are < 0.019, < 0.014, and -0.020, respectively. 
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4.3. How did the cores fragment? 

Recent numerical simulations and observations support the hypothesis that the fragmen- 
tation of molecular cloud cores is the main mechanism for the formation of binary /multiple 
stellar systems, although the exact when, where, why, and how are still under debate (see 
reviews by Bodenheimer et al. 2000, Tohline 2002, and Goodwin et al. 2007). In this section, 
we try to examine the origin of the sub-cores in both CG30 and BHR71, i.e., whether they 
formed by initial cloud fragmentation prior to protostellar collapse or by prompt rotational 
fragmentation of a single core after the initial collapse. 

In CG30, our previous single-dish submm maps have shown a large-scale hourglass- 
shaped common envelope around the two sub-cores (Henning et al. 2001; see Fig. la). The 
separation between the sub-cores is ~ 8700 AU, which is roughly two times the typical Jeans 

length [i?jeans = 0.19 

P*"(laff)^ ( io4cm-3 ) ^' Stahler & Palla 2004] in prestellar cores 
(~4000AU at T = 10 K and riHa = lO^cm"^). The radial velocity difference between the 
two sub-cores is ~0.16kms~^ (see Table 3). If we assume that the total binary mass is 
1.4 X 1.36 Mq (see Table 2; the factor 1.36 accounting for He and heavier elements) and the 
orbit is perpendicular to the plane of sky, the orbit velocity difference in a bound system 
with the separation of 8700 AU should be ~0.44kms~^, about three times larger than the 
observed value. Furthermore, from this observed velocity difference, we estimate the /^rot of 
~ 0.008 for the large-scale cloud core which contains the two sub-cores (radius ~8000AU; 
see Fig. 2a). This Prot is less than the typical boundary suggested by a series of numerical 
simulations (see e.g.. Boss 1999 and Machida et al. 2005) for rotational fragmentation. Based 
on the morphology and velocity structure, we suggest that the two sub-cores in CG 30 were 
formed by initial fragment at ioE(^ of a large-scale filamentary prestellar core. 

In BHR71, the two sub-cores have a separation of ~3400AU (less than the typical 
Jeans length) and are also surrounded by a large common envelope (Fig. lb). Unfortunately, 
the observed velocity structure is mainly associated with IRSl and kinematic information of 
IRS2 is missing. Here we can only speculate on the basis of separation that the two sub-cores 
could be formed by prompt rotational fragmentation of a collapsing protostellar core. 



The basic idea of this initial fragmentation is that the collapse is initiated in a large-scale molecular cloud 
core which contains multiple Jeans masses in a weakly condensed configuration, e.g., a prolate or filamentary 
Gaussian distribution with several Jeans masses along the long axis and one Jeans mass across the short axis; 
with some initial angular momentum, provided by either slow rotation (Bonnell et al. 1991) or turbulence 
(Goodwin et al. 2007), the large cloud core fragments at Jeans scale into several dense cores, in which the 
separate protostellar collapse then starts and proceeds more quickly than across the whole structure (see 
e.g., Mundy et al. 2001 and reference therein). 
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Numerical simulations also predict that the material collapses along the magnetic field 
lines while the fragmentation occurs in a plane perpendicular to the magnetic field. This is 
supported by our previous submm polarimetric observations towards CG 30 (Henning et al. 
2001). In contrast to a simple assumption that the angular momenta of two components will 
be parallel in the fragmentation, we find that the outfiows, and hence the angular momenta 
(assumed to be parallel to the outfiows), of the sub-cores are not aligned, neither in CG30, 
nor in BHR71. This phenomenon is also found in other binary protostars studied recently, 
like e.g., CB230 and L723 (Launhardt 2004; Launhardt et al., in prep.). This could mean 
that during core fragmentation the initial angular momentum is not evenly (in value and 
direction) divided between the sub-cores, although the mean direction is preserved all the 
time. 



4.4. N2H+ vs. Dust vs. CO 

Prom our observations towards Class protostars conducted at OVRO (Paper I), ATCA 
(this work), and IRAM-PdBI (Chen et al., in prep.), we find that in most objects the mm 
continuum source lies within the half maximum level of the N2H"'" emission. This good 
general agreement indicates that N2H''" is spatially associated with thermal dust in dense 
protostellar cores and cannot be significantly depleted like, e.g., CO and CS (see Bergin et 
al. 2001 and Caselh et al. 1999). 

Pig. 10 shows that the dust mass (converted into hydrogen gas mass) is in general cor- 
related with both the N2H"'" gas mass and the virial mass (both derived from the N2H"'" 
emission). However, there is a significant scatter in both correlations, indicating that the 
agreement holds only within a factor of 2 to 2.5 (1 a scatter). This could be due to the fact 
that the mm dust continuum emission traces mainly the dense structures (e.g., inner enve- 
lope or disk), while N2H+ emission traces the larger-scale envelope (see e.g.. Pigs. la&2a), 
and hence refiects different masses on different scales. We also note that the N2H+ gas mass 
depends on the specific source morphology and chemistry since it is quickly destroyed where 
CO is released from dust grain into the gas phase (see below). The estimated virial mass 
has also significant uncertainties because several sources are driving bipolar outfiows and are 
probably no longer in virial equilibrium. 

On the other hand, we also find that in most objects the morphology of the N2H+ emis- 
sion is directly related to the jet/outfiow actions. Por example, in BHR71, two N2H"'" cores, 
located to the east and west of the outfiow-driving source IRSl, are rotating perpendicular 
to the outfiow axis, and there is no N2H+ emission detected at the origin and along the 
large-scale CO outfiow. These features suggest a large N2H+ hole has been formed and the 
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two cores may be the remnant of a N2H"'" envelopq^ It is likely that a large amount of 
N2H+ in the way of the outflow has been depleted by CO molecules, which is one of the 
main destroyers of N2H"'" in the gas phase (Aikawa et al. 2001). For this reason, we think 
that the emission at the position of IRS2 (see Fig. 2b) is part of the structure around IRSl 
and does not originate from IRS2. 

Based on the observational results, we speculate that there are three stages of the 
interaction between N2H+ and jets/outfiows. (1) When jets are ejected from a protostar, 
N2H+ molecules in the envelope are entrained and show a jet-like morphology in the images, 
like, e.g., L723 VLA2 (Paper I) and CG30N (this work). (2) Molecular outflows, following 
the jets, release CO from grain surfaces back into the gas phase and start destroying the 
N2H+ molecules on the way, leading to the observed hourglass structure perpendicular to the 
CO outflow axis, like e.g., IRAS 03282+3035, IRAS 04166+2706, and CB224 (see Paper I). 
(3) Large N2H"'" holes form in the envelopes, like e.g., in BHR71 IRSl (this work), L483 
(J0rgensen 2004), and IRAM 04191 (Belloche & Andre 2004). However, there seems to be 
no clear correlation between this N2H"'"/jet scenario and standard evolutionary scenario from 
Class to Class I. For example, IRAM 04191 is a young Class protostar but appears in the 
last stage, while CG30N is a Class I object but appears in the first stage. We speculate that 
the appearance of the N2H"'" emission is strongly affected by outflow-envelope interaction, 
which depends on the specific envelope morphology and source multiplicity properties. 



5. SUMMARY 

We have presented ATCA and Spitzer observations of the two isolated protostellar 
double cores CG30 and BHR71 in the southern sky. The main results of this work are 
summarized as follows: 

(1) The 3 mm dust continuum image of CG30 resolves two compact sources with a 
separation of ~ 21" (8400 AU). In BHR71, one strong dust continuum source is detected 
at the position of mid-infrared source IRSl, while only weak emission is detected from 
the secondary mid-infrared source IRS2. The separation between IRSl and IRS2 is ~ 17" 
(3400 AU). Assuming optically thin dust emission, we derive hydrogen gas masses of 1.1 Mq 
and 0.33 Mq for northern and southern sources in CG30, and 2.1 Mq and 0.05 Mq for IRSl 
and IRS2 sources in BHR71. 



^'^For similar cases see low-mass protostars L483 (J0rgensen 2004) and IRAM 04191 (Belloche & Andre 
2004). 
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(2) N2H+ (1-0) emission is detected in both CG30 and BHR71. In CG30, the two 
dust continuum sources are directly associated with N2H"'" cores. In BHR71, two N2H"'" cores 
are around the primary dust continuum source, probably part of one large envelope, but no 
N2H+ is detected at the position of the dust source. The secondary IR source is not detected 
in N2H+. 

(3) The excitation temperatures of the N2H+ line are 4.7 — 6.8 K for CG 30 and 3.9 — 4.4 K 
for BIIR71. The FWHM radii of N2H^ cores range from 730 to 1700 AU. The average frac- 
tional abundances of N2H+, derived from the ratio of N2H"'" gas mass to virial mass, is 
~3.0 X 10"^'', which is consistent with the results obtained in our previous study of the 
cores in northern sky. The observed mean N2H"'" line widths are ~0.5kms^^ for CG30 and 
~0.3kms~^ for BHR71. The line widths are roughly constant within the interiors of the 
cores and large hne widths only occur at the edges of the cores. The derived virial masses 
of the N2H+ cores range from 0.1 to 0.6 M©. 

(4) We derive the N2H+ radial velocity fields for GG 30 and BHR 71. The two N2H+ cores 
in GG 30 show systematic velocity gradients of ~ 24.4 km s~^ pc~^ and ~ 17.8 km pc~^ that 
are parallel to the outfiow directions and could be affected by the outfiows. In BHR 71, a 
systematic velocity gradient of ~ 7.8kms~^ pc~^ across the two cores is perpendicular to the 
large-scale outfiow and could be explained by rotation. 

(5) Assuming that the observed velocity gradients are due to core rotation (if perpen- 
dicular to outfiow) or place an upper limit to rotation (if parallel to outfiow), we estimate 
specific angular momenta of < 0.30, < 0.35, and ~ 0.51 x 10"^ kms"^ pc for GG 30N, GG 30S, 

and BHR 71, respectively. The ratios for the rotational energy to the gravitational potential 
energy for CG30N, GG30S, and BHR 71 are estimated to be < 0.019, < 0.014, and -0.020, 
respectively. 

(6) Infrared emission from both sub-cores in both CG30 and BHR 71 is detected at 
Spitzer IRAC bands and MIPS bands. Each source is driving its own outflow, as seen in 
the shock-excited 4.5 fim infrared images. CG 30N is associated with a Herbig-Haro flow, 
while the southern source is driving a large bipolar jet. In BHR 71, both IRSl and IRS2 are 
associated with Herbig-Haro objects and driving bipolar jets which coincide spatially with 
the GO outfiows. 

(7) By fitting the spectral energy distributions, we derive the dust temperature, bolo- 
metric temperature, and bolometric luminosity of the sources. We flnd that CG 30N is a 
Class I object while the southern source is a Class protostar. In BHR 71, the properties of 
IRSl resemble a Glass O/I transition object, while IRS2 is a Glass protostar. We speculate 
that the sources may nevertheless be coeval but that this evolutionary discrepancy is due to 
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the effects of ralative inclinations. 

(8) Based on the morphologies and velocity structures, we suggest that the double cores 
in CG 30 were formed by initial fragmentation of a filamentary prestellar core, while BHR 71 
may originate from rotational fragmentation of a single collapsing protostellar core. We also 
find that the angular momenta of the sub-cores are not aligned in either pair of sources. 

(9) Our observations conducted at OVRO and ATCA show a close correlation between 
thermal dust emission and N2H+. The N2H+ emission in most sources is spatially associated 
and quantitatively correlated with the dust continuum emission. However, we also find 
a strong relationship between the morphology of the N2H+ emission and the jet/outfiow 
actions. Outflows first seem to entrain N2H"'' and then gradually destroy it, which leads to 
the observed jet-like, hourglass-shaped intensity maps and N2H"'" hole. 

We thank the anonymous referee for many helpful comments and suggestions. The 
Australia Telescope Compact Array is part of the Australia Telescope, which is funded by 
the Commonwealth of Australia for operation as a national facility managed by CSIRO. 
Wc thank the ATCA staff for technical support during the observations. We also thank 
A. Goodman for fruitful discussions and providing the VFIT routine. 
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Table 1. Target list and summary of observations 



Object 


IRAS 


R.A. & Dec. (J2000)'' Distance Array 


HPBW'' 


rms"^ 




Name 


Source 


[h : m : s, ° : ' : "] [pc] configuration 


[arcsec] [mjy/beam] 


CG30 


08076-3556 


08:09:33.0, -36:05:01.00 400 H75+H168 


4.6x3.3/4.6x3.2 


65/0.5 




BHR71 


11590-6452 


12:01:36.5, -65:08:49.49 200 H75+H168 


3.6x2.9/3.9x3.1 


on /o n 




'^Reference position for figures and tables in the paper (except Spitzer images) 








''Synthesized FWHM beam sizes at N2H+(1 — 0) hne / 3 mm dust continuum with robust weighting 1. 






''la noises at N2H+(1- 


0) line / 3 mm dust continuum. 










Table 2. 


3 mm dust continuum results for CG 30 and BHR 71 






Source 


R.A." 


Dec." FWHM sizes" 


Mh (nn)'' 








(J2000) 


(J2000) [mjy] maj.xmin. P. A. 


[Mq] [xlO'^cm- 


-3] [X 


1023cm-2l 


CG30N 


08:09:33.12 


-36:04:58.12 15.8±3.2 5'.'lx3'.'l 89±7° 


1.10±0.26 1.11 




4.51 


CG30S 


08:09:32.67 


-36:05:19.09 6.0±1.3 4y8x3^'l 74±15° 


0.33±0.10 0.37 




1.45 


BHR71 IRSl 


12:01:36.81 


-65:08:49.22 140±28 7^'8x7^'l -73±20° 


2.12±0.41 2.64 




9.94 


BHR71 IRS2 


12:01:34.09 


-65:08:47.36 2.8±2.1 2^'6x2^'l 76±40° 


0.05±0.02 2.18 




2.58 


^Center position and FWHM 


sizes of the continuum sources derived from Gaussian 


uv plane fitting. 






''Assuming a 


spherical morphology for the objects, the mean volume density of hydrogen atoms nn = 


n(H) + 


2n(H2) was 


calculated by riH 


= Mu/mnV, with V ~ n/6{6sD)^ being volume. 









'^The hydrogen column density ATh = Ar(H) + 2 Ar(H2) was derived from the flux densities by A^h = 
where Qg is the solid angle of the objects and mn is the proton mass. 



Table 3. Observing parameters from N2H+ (1 — 0) spectra fitting 



Source 


[km s-i] 


Ai;" 
[km s-l] 


not" 


T' 0, 
ex 

[K] 


All ^ 
[km s-l] 


[km s-i] 


CG30N 


6.64±0.02 


0.53±0.03 


l.OiO.l 


4.66±0.09 


0.51±0.01 


0.47 


CG30S 


6.48±0.01 


0.52±0.02 


1.4±0.1 


6.81±0.05 


0.52±0.01 


0.48 


BHR 71 E 


-4.35±0.02 


0.44±0.10 


1.9±0.4 


4.44±0.13 


0.28±0.01 


0.20 


BHR 71 W 


-4.42±0.02 


0.38±0.06 


2.5±0.2 


3.91±0.05 


0.33±0.01 


0.27 



''Value at the intensity peak. The error represents 1 cr error in the hyperfine fitting. 

''Mean line width obtained through Gaussian fitting to the distribution of line widths 

versus solid angle areas. 

"^Non-thermal line width at the given dust temperature (see Table 7). 
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Table 4. Volume size, density, and mass of N2H+ cores 



Source 


R 

[AU] 


Mvir 


[xlO« cm-3] 


JV(N2H+) 
[XIOI2 cm-2] 


[x 10-10 Mq] 


X(N2H+) 
[x 10-10] 


CG30N 


1300 


0.38 


6.2 


1.29 


0.51 


1.83 


CG30S 


1650 


0.55 


4.5 


3.42 


1.94 


4.85 


BHR71E 


960 


0.13 


5.3 


1.95 


0.34 


3.63 


BHR71W 


730 


0.11 


10.3 


1.87 


0.21 


2.52 


BHR71 IRSl"" 


3000 


0.49 


0.7 


1.90 


3.07 


8.55 



^Assuming that the two N2H+ cores found in BHR71 arc part of a physical structure around 
IRSl with a radius of ^SOOOAU (see Fig. 2b) and a mean Une width of 0.33kms— ^ (see Table 3). 



Table 5. Velocity gradients and specific angular momentum 





mean velocity 


9 


©9 




J/M 




Source 


[km s— 1] 


[km s^l pc^i] 


[degree] 


[km s^i] 


[x 10-^ km s^i pc] 




CG30N 


6.62 


<24.4±0.2 


-79.7±0.4 


0.32 


< 0.30 


<0.019 


CG30S 


6.45 


<17.8±0.2 


33.9±0.3 


0.29 


< 0.35 


<0.014 


BHR71 IRSl'' 


-4.39 


~7.8±0.5 


-104±2.0 


0.23 


~ 0.51 


~0.020 



*East of north in the direction of increasing velocity 
''The same assumption as in Table 4. 



Table 6. Spitzer flux densities of CG 30 and BHR71" 





R.A.'' 


Dec.'' 


5(3.6 fj,m) 


5(4.5 ^lm) 


5(5.8 /im) 


5(8.0 /im) 


5(24 /im) 


5(70 /im) 


Source 


(J2000) 


(J2000) 


[mjy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


CG30N 


08:09:33.20 


-36:04:58.17 


55.7±1.4 


123.4±2.1 


256.1±3.0 


395.8±3.8 


3400±100 


8700±430 


CG30S 


08:09:32.68 


-36:05:20.38 


6.7±0.5 


16.8±0.8 


19.6±0.8 


9.6±0.6 


50±5 


4200±340 


BHR71 IRSl 


12:01:36.57 


-65:08:49.52 


32.4±1.1 


82.4±1.7 


123.3±2.1 


210.2±2.8 


5000±300 


84000±800 


BHR71 IRS2 


12:01:34.05 


-65:08:47.03 


4.5±0.4 


12.4±0.7 


15.4±0.7 


9.3±0.6 


90±30 





^Flux densities in the IRAC and MIPS bands were measured using IRAF APPHOT and GILDAS, respectively (see §2.2). 
''Peak position of infrared sources measured at the IRAC band 3 (5.8 //m). 



Table 7. Fitting results of the spectral energy distribution 



Source 


^dust 






-^submm 


-^submm /^bol 


Classification 




[K] 


[K] 


[Lq] 




[%] 




CG30N 


22 


102 


13.6±0.8 


0.49±0.10 


3.6 


Class I 


CG30S 


27 


37 


4.3±0.5 


0.32±0.05 
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Fig. 1. — (a) 3 mm dust continuum image of CG30. Contours start at ~ 3 cr (Icr^O.SmJy) 
with steps of ~ 2 o". The grey dashed contour represents the half-maximum level of the 
850 /im emission observed with SCUBA (Henning et al. 2001). (b) 3mm dust continuum 
image of BHR 71. Contours start at~3o'(lo"~2 mJy) with steps of ~ 2 o". Crosses mark the 
positions of the Spitzer MIR sources. The grey dashed contour represents the half-maximum 
level of the 1.2 mm emission observed with SEST (B97). Synthesized ATCA beams are shown 
as grey ovals in bottom right corners. 
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Fig. 2. — (a) Image of the N2H+ (1 — 0) intensity integrated over the seven hyperfine compo- 
nents for CG 30. Contours start at ~ 3 cr (1 cr ~ 60 mJy) with steps of ~ 2 cr. The arrows show 
the directions of protosteUar jets (see §3.3). The grey dashed contours represent the 50%, 
75%, and 99% levels of the N2H+ (1 — 0) emission observed with Mopra single-dish telescope 
(P. Barnes et al. in prep.), (b) The same for BHR71 (l(T~20mJy). The solid and dashed 
arrows show the directions of the blue-shifted and red-shifted CO outflows (see §3.3). The 
crosses in both images represent the peaks of 3 mm dust continuum emission. Synthesized 
ATCA beams are shown as grey ovals. 
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Fig. 3. — N2H"'" spectra at the peak positions of the two cores in CG30 (left) and BHR71 
(right). Thin dotted curves show the results of hyperfine structure line fitting. Fit parameters 
are given in Table 3. [See the electronic edition of the Journal for a color version of this 
figure] 
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Fig. 4. — N2H"'" velocity field maps of CG30 (left) and BHR71 (right). Contours are same 
as in Fig. 2. The white arrows in CG SON and CG 308 show the directions of the protostellar 
jets. The red and blue arrows in the BHR71 map show the directions of CO outflows, while 
the white arrow shows the direction of the gradient across the two main N2H"'' clumps. 
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Fig. 5.— Spatial distribution of N2H+ line widths in CG30 (left) and BHR71 (right), as 
derived from the HFS line fitting. Contours and symbols are the same as in Fig. 2. 
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Fig. 6. — Distribution of N2H"'" line widths versus solid angle areas for CG30 (top) and 
BHR71 (bottom). Black solid curves and numbers show the results of Gaussian fitting to 
the distributions. 
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Fig. 7 — Spitzer images of CG30. (a) Spitzer IRAC band 2 (4.5 /im) image of CG30. Sources CG30N 
and CG SOS are labeled as "N" and "S" , respectively. Dashed arrows show the directions of the protostellar 
jets; (b) IRAC band 2 image overlaid with the ATCA 3 mm dust continuum contours (reference position 
at R.A.=08:09:04.082, DEC=-36:00:53.53, J2000); (c) Same, but overlaid with the ATCA N2H+ intensity 
contours; (d) Spitzer MIPS 1 (24 /im) image of CG30, overlaid with the ATCA 3 mm dust continuum 
contours; (e) Spitzer MIPS 2 (70 /xm) image of CG 30, overlaid with the ATCA 3 mm dust continuum 
contours. 
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Fig. 8. — The same as Fig. 7, but for BHR71 (reference position at R.A. =12:01:36. 349, 
DEC=-65:08:37.50, J2000). 
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Fig. 9. — Spectral energy distribution of CG30N (up left), CG30S (up right), BHR71 IRSl (bottom left), 
and BHR71 IRS2 (bottom right). Error bars (Icr) are indicated for all data points, but are mostly smaller 
than the symbol sizes. Open squares represent IRAS data points, where flux densities are divided into two 
sub-cores with ratios assumed in §4.1. While most data points represent total fluxes, the 3 mm fluxes were 
measured from interferometric maps which resolved out the envelope and thus represent lower limits only. 
Solid lines show the best-fit for all points at A > 100 /im using a grey-body model. Dashed lines at A < 100 /im 
show the simple logarithmic interpolation used to derive the luminosity. The fitting results are summarized 
in Table 7. 
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Fig. 10. — (a) N2H+ gas mass (derived from N2H+ line emission) versus hydrogen gas mass (derived from 
3 mm dust continuum emission; the dust continuum data of OVRO sample are taken from Launhardt et al. 
in prep.), and (b) virial mass (derived from N2H+ line emission) versus hydrogen gas mass for protostellar 
cores studied in Paper I and this work. 



